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1.  INTRODUCTION 


As  indicated  in  Figure  1,  the  problem  we  are  addressing  is  that  of  battlefield 
clutter  due  to  the  presence  of  aerosols  and  atmospheric  motion  produced  by  turbulence 
and  rapidly  moving  hot  air.  We  have  a  variety  of  aerosol  sources  such  as  dirt  kicked 
up  by  explosions,  dust  clouds  due  to  wind  storms  and  vehicular  motion,  haze  and  fog, 
smoke  from  fires,  and  screening  smokes.  Each  of  these  aerosol  sources  contributes 
serious  degradation  to  the  operation  of  electro-optical  systems.  (Xir  work  has  been 
directed  towards  describing  and  developing  an  understanding  of  the  various  aerosol 
source  terms  and  their  subsequent  diffusion  in  the  atmosphere.  From  this,  we  have 
been  developing  a  disturbed  IR  transmission  (DIRTRAN)  model  to  describe  the  degra¬ 
dation  effects  produced  by  battlefield  clutter.^  **)  Validation  of  this  model  is  being 
performed  with  the  aid  of  data  obtained  at  several  recent  Army  infrared  obscuration 
field  trials. 

The  goal  of  the  Aerodyne  DIRTRAN  model  is  two-fold:  (1)  predict  the  history  of 
aerosol  clouds  produced  in  a  battlefield  environment  (specifically  for  this  report,  we 
discuss  dust  clouds  produced  by  explosions  in  soil  in  different  geographic  scenarios); 
(2)  predict  the  time-dependent  obscuration  at  several  wavelengths,  from  the  visible  to 
the  near  IR  to  mm  wavelengths.  In  the  following  three  sections,  we  discuss  three 
major  elements  of  the  DIRTRAN  model,  as  schematically  shown  in  Figure  2.  There 
is,  first,  the  aerosol  source  term,  hi  the  case  of  artillery  explosions,  there  are 
bomb  crater  dynamics  and  soil  characteristics  which  must  be  taken  into  account  in 
order  to  have  an  understanding  of  the  particle  ejecta.  Next,  to  adequately  account  for 
atmospheric  diffusion,  the  effects  of  wind,  turbulence,  and  gravity  must  be  included. 
Finally,  the  obscuration  clutter  produced  in  an  EO  detection  system  will  be  the  result 
of  variable  extinction,  scattering,  absorption,  and  thermal  emission  from  the  time- 
dependent  aerosol  clouds. 


BATTLEFIELD  CLUTTER 


r  .v».«p<! 


Laser  Designator 
Target 


Figure  1.  Example  of  a  Clutter  in  a  Battlefield  Scenario.  This  results 
from  the  presence  of  aerosols  and  turbulence  from  a  variety 
of  sources;  dust  from  explosions,  vehicles,  and  dust  storms; 
screening  smokes;  fire  smoke;  fog  and  haze. 
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ELEMENTS  OF  AERODYNE  DISTURBED  INFRARED  TRANSMISSION  MODEL  -  DIRTRAN 


Source  Atmospheric  i 

Term  - -  - - Diffusion  — -I- - Clut1er 


Bomb  Crater 
Dynamics 


Soil 

Characteristics 


Particle 

Ejecta 


Wind,  Lapse  Rate, 
(Turbulence) 
and  Gravity 


Differential  Particle 
,  Entrainment,  Diffu¬ 
sion  and  Settling 
(Dust  Cloud  Dispersal) 


>  Obscuration 


Optical 

Transmission 


Figure  2.  Elements  of  the  Aerodyne  DIRTRAN  (Disturbed  Infrared  Trans¬ 
mission)  Model.  The  aerosol  source  term  (e.  g. ,  dust  resulting 
from  an  explosion  in  soil)  is  first  described.  The  aerosols  are 
then  subjected  to  atmospheric  diffusion  resulting  from  wind, 
turbulence,  and  gravitational  settling.  Finally,  the  clutter 
(obscuration  and  scintillation)  of  optical  radiation  is  computed. 
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2.  AEROSOL  SOURCE  TERM;  EXPLOSION- PRODUCED  DUST 


The  amount  of  dust  generated  by  an  explosion  in  soil  is  obtained  from  the  size 
of  the  crater  predicted  in  accord  with  Ref.  (7)  as  shown  in  Figures  3  and  4.  The 
particle  size  distribution  function  in  the  ejecta  for  the  smaller  radii  (<100  fim)  is 
assumed  to  be  log-normal  with  a  mean  radius,  standard  deviation,  and  mass  loading 
which  will  be  determined  from  a  systematic  analysis  of  data.  It  is  noted  that  "the 
process  of  comminution  of  the  crater  material  has  not  been  described  either  theoret¬ 
ically  or  experimentally"  (Ref.  (8) ).  A  comparison  of  soil  analysis  and  dust  cloud 
particulate  data  will  provide  insights  into  this  problem.  Comparable  uncertainty  also 
exists  with  regard  to  the  initial  velocities  of  the  ejecta  (Refs.  (8)  and  (9) ).  Aero¬ 
dynamic  deceleration  causes  the  ejecta  to  come  to  rest  relative  to  the  air  and  thus 
determines  the  initial  conditions  for  atmospheric  diffusion  of  the  dust.  This  aspect 
of  the  problem  is  important  in  determining  the  partitioning  of  soil  ejecta  inside  and 
outside  the  hot  gases  as  discussed  below  in  Section  2C.  These  will  be  treated  paramet¬ 
rically  and  compared  with  data. 


Effect  of  Vegetative  Cover  on  Production  of  Dust  Clouds 


Figure  5  shows  the  results  of  transmission  tests  in  the  8-12  ^m  band  performed 
by  the  Army  Night  Vision  &  Electro- Optics  Laboratory  (NV&EOL),  Ft.  Belvoir,  Va. 

This  test  —  known  as  TAPATS  (Threat  Artillery  Preparation  Against  Tank  Sights) 
was  conducted  in  July  1978  at  Ft.  Knox,  Ky.  The  transmission  data  shown  are  the 
results  of  a  tank  night  sight  viewing  through  a  155  mm  HE  barrage  of  intensity 
0. 8  rounds /sec /km.  Notable  in  the  figure  is  the  fact  that  from  the  first- round  impact 
to  several  tens  of  seconds  after  the  barrage  ended,  the  transmission  was  essentially 
zero.  That  is,  there  were  no  visibility  holes  evident  for  the  entire  length  of  the  bar¬ 
rage.  In  Figure  C  is  plotted  the  transmission  during  a  barrage  conducted  at  Grafenwohr, 
Germany  tn  November  1978  by  NV&EOL.  In  this  case,  the  barrage  had  an  intensity  of 
1.  5  rounds /sec /km.  Here,  attenuation  due  to  background  fog  and  natural  aerosols 
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Figure  4.  Nomograph  Similar  to  Figure  3  But  Used  for  Determining 
the  Apparent  Crater  Depth  d  in  Various  Media  (Ref.  7). 


(Ft,  Knox,  Ky.) 


Wavelength  Band:  8  -  12  hm 
Artillery:  155  mm  HE 
Barrage:  0.8  Rounds/ km/slc 
(no  visibility  holes  evident)  Impact  Area:  100  m  x  300  m 


Atmosphere:  Stable 

WlNDSPEED:  1  KNOT 


1 1 me  After  Beginning  of  Barrage  (sec) 


Figure  5.  Transmission  in  the  8- 12  Band  During  a  155  mm  HE  Barrage  at 
Ft.  Knox,  Ky.  Far  this  barrage  intensity  and  atmospheric  condition, 
essentially  zero  transmission  was  observed  for  the  duration  of  the 
barrage  and  several  tens  of  seconds  after  it.  (Data  from  the  TAPATS 
test  conducted  by  the  Night  Vision  &  Electro- Optics  Laboratory, 

Ft.  Belvoir,  Va. ) 
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Transmission  (Z) 


6ML.11 


0  30  60  90  120  150  180  210  290 


Time  After  Beginning  of  Barrage  (sec) 


Figure  6.  Transmission  in  the  8-12  ^m  and  0.  8-  1. 1  ^m  Band  During  a  155  mm  HE 
Barrage  at  Grafenwohr,  Germany.  In  this  situation,  visibility  holes  in 
excess  of  20-40%  transmission  occurred  at  numerous  times  during  the 
barrage.  (Data  from  the  Graf  II  test  conducted  by  the  Night  Vision  & 
Electro-Optics  Laboratory,  Ft.  Belvoir,  Va. ) 
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(as  determined  prior  to  the  barrage)  have  been  removed  from  the  transmission  data, 
so  that  what  is  shown  in  the  figure  is  strictly  the  changes  in  transmission  due  to  the 
barrage.  It  is  evident  from  the  figure  that  the  transmission  was  considerably  in 
excess  of  20  to  40%  for  several  tens  of  seconds  at  several  different  times  during  the 
barrage.  That  is,  even  with  a  barrage  intensity  nearly  twice  as  much  as  the  TAPATS 
test,  the  Graf  II  test  transmission  was  considerably  better. 

In  order  to  ascertain  the  reasons  for  the  differences  in  transmission  between 
the  Graf  n  and  TAPATS  tests,  we  looked  at  a  variety  of  possible  circumstances.  One 
of  the  first  things  examined  was  the  effect  of  soil  moisture.  Table  I  is  a  summary'  of 
the  moisture  content  of  the  soil  for  the  TAPATS  and  Graf  II  tests.  It  may  be  noted 
that  the  average  moisture  content  for  Zone  4,  Grafenwohr  (corresponding  to  Figure  6) 
was  approximately  20%;  the  moisture  content  for  the  TAPATS  test  was  on  the  order 
of  16%.  That  is,  the  moisture  content  for  both  tests  were  comparable.  We  were 
thus  led  to  conclude  that  soil  moisture  content  was  not  the  primary  mechanism  for 
describing  the  difference  in  transmission  between  the  TAPATS  and  Graf  II  tests. 


TABLE  I 

SOIL  MOISTURE  CONTENT 


Test 

Location 

Date 

Moisture  Content 

TAPATS 

Ft.  Knox,  Ky. 

July  1978 

16% 

Graf  II 

Grafenwohr,  Germ 

Nov.  1978 

13  -  27%  (Zone  4) 

6.  5  -  15%  (Zone  1) 

We  next  looked  at  photographs  showing  the  size  of  the  craters  produced  in  the 
Graf  II  tests.  Immediately  evident  is  the  fact  that  the  depth  of  the  grassy  sod  is  on 
the  order  of  25  cm,  as  shown  in  Figure  7.  Also  evident  in  the  figure  is  the  existence 
of  large  clumps  of  sod  which  have  been  simply  rolled  over  due  to  the  explosion  of  the 
155  mm  high  explosive  round.  Figure  8  is  a  diagram  showing  the  effect  of  sod 
on  crater  volume.  It  is  drawn  approximately  according  to  scale  for  the  crater  depth 
produced  by  a  155  mm  HE  round.  It  may  be  seen  that  the  sod  takes  up  the  upper 


Figure  7.  Graf  II  Test  155  mm  HE  Hound  Craters.  Depth  of  vegetative 
sod  cover  is  approximately  25  cm,  comprising  a  large  portion 
of  the  crater  volume. 


Figure  8.  Effect  of  Sod  (Vegetative  Cover)  on  Crater  Volume.  The  sod 
depth  is  drawn  here  roughly  to  scale  for  the  Graf  II  test  for 
the  155  mm  HE  round  craters.  The  volume  of  soil  (which  is 
available  for  production  of  entrained  dust)  is  seen  to  be  only 
a  small  fraction  of  the  total  crater  volume. 


portion  of  the  crater  (where  the  greatest  crater  volume  occurs).  Only  a  relatively 
small  portion  of  the  total  crater  volume  contains  soil  particles  which  might  be  entrained 
in  the  atmosphere.  To  see  this  effect  mathematically,  if  we  assume  the  crater  volume 
is  conical  in  shape,  then  that  part  of  the  crater  volume  Vc  which  contains  only  soil 
(i.  e. ,  no  sod)  is  given  by  the  fallowing  equation: 


vc  =  i  [t )  <d~  -  <u 


where  rc  and  dc  are  the  crater  radius  and  diameter  respectively  (as  determined  in 
accord  with  Ref.  (7) )  and  where  d  is  the  sod  depth.  It  may  be  noted  here  that  the 
effect  of  sod  depth  comes  in  as  the  cube  of  the  difference  between  the  crater  depth  and 
the  sod  depth,  i.  e. ,  the  soil-only  crater  volume  is  a  very  sensitive  function  of  sod 
depth. 
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The  Impact  zone  for  the  Ft.  Knox  TAPATS  test  is  being  continually  churned  up 
by  test  barrages.  That  is,  there  is  no  substantial  vegetative  or  sod  cover  because  it 
is  continually  being  depleted  by  the  explosion  of  HE  rounds.  Thus,  we  are  led  to  the 
conclusion  that  the  differences  in  the  obscuration  between  the  TAPATS  and  Graf  II 
data  may  be  attributed  to  the  presence  of  vegetative  soil  cover  (sod)  in  the  Germany 
test  site,  a  sod  thickness  sufficient  to  prevent  much  dust  entrainment,  allowing  a 
considerable  number  of  visibility  holes  to  occur  in  the  transmission  data. 

2 C.  partitioning  of  Particles  Inside  and  Outside  of  the  Explosion  Fireball 

Figure  9  is  a  photograph  which  was  taken  during  the  Dusty  Infrared  Test  I 
(DIRT  I)  performed  by  the  Atmospheric  Sciences  Laboratory,  White  Sands  Missile 
Range  in  October  1978.  Evident  in  the  photograph  is  the  buoyantly  rising  cone-shaped 
or  v-shaped  hot  ball  of  gases  and  dust  produced  by  the  explosion,  i.  e. ,  the  hot  gases 
and  dust  first  rise  into  the  atmosphere  before  being  dominated  by  wind  diffusion.  Sur¬ 
rounding  this  hot  gas  at  the  base  is  a  skirt  of  dust  that  is  not  buoyantly  rising  and  which 
is  immediately  subjected  to  wind  diffusion.  Thus,  both  aspects  of  this  aerosol  source 
term  must  be  taken  into  account  with  regard  to  atmospheric  diffusion.  Figure  10  is 
also  taken  from  the  DIRT  I  test  series.  In  this  case,  the  explosion  has  occurred 
several  seconds  prior  to  the  taking  of  the  photograph.  It  may  be  noted  that  the  hot 
fireball  gases  and  dust  have  risen  to  a  higher  altitude  and  that  there  is  a  puff  at  the 
top  of  the  cloud.  Noteworthy  also  is  the  fact  that  the  skirt  of  dust  at  the  base  of  the 
cloud  is  still  present.  This  will  have  a  significant  effect  on  ground-based  EO  systems. 

2D.  Initial  Hot  Gas  Cloud 

The  formation  and  dispersion  of  dust  clouds  produced  by  artillery  explosions 
involve  two  distinct  sequential  phases  characterized  by  two  drastically  different  time 
scales.  In  the  first  phase,  formation  of  the  crater  and  the  pressure/velocity  equili¬ 
bration  of  the  detonation  products  as  well  as  the  ejecta  occur  within  times  on  the 
-  3  -2 

order  of  10  to  10  sec.  The  DIRTRAN  model  relies  on  existing  methodology  for 
the  state  and  size  of  the  cloud  at  the  end  of  this  phase.  Specifically,  the  size  and 
temperature  of  the  pressure- equilibrated  hot  gases  are  determined  from  well 
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established  correlations  of  hydrodynamic  blast  calculations  (Ret  (10)).  The  second 

involves  subsequent  buoyancy  and  diffusion-controlled  dispersion  which  takes  place  on 

1  2 

time  scales  on  the  order  of  10  to  10  sec  and  is  described  in  Section  3. 


Figure  9.  Example  of  Entrainment  partitioning  of  Explosion- produced  Soil 
Ejecta  Inside  and  Outside  of  Hot  Gases  From  Explosion. 

(Data  from  DIRT  I  Test  series  conducted  at  White  Sands  Missile 
Range  in  October  1978  by  the  Army  Atmospheric  Sciences 
Laboratory. ) 


Figure  10.  Example  of  Buoyant  Rise  Effects  on  the  Dispersion  of  Explosion- 
Produced  Dust  Clouds 

Also  evident  is  the  effect  of  a  mild  wind  which  causes  the  rising 
cloud  to  follow  a  slanted  path.  In  addition,  dust  which  was  not 
entrained  in  the  hot  gases  is  still  present  near  the  explosion 
point.  (Data  from  the  DIRT  I  Test  series.) 


3.  ATMOSPHERIC  DIFFUSION 


3A.  Buoyant  Cloud  Rise  and  Diffusion 

The  model  presented  in  Ref.  (11)  for  instantaneous  point  releases  are  used. 
Modifications  are  introduced  to  accommodate  the  finite  initial  cloud  size  and  arbitrary 
vertical  temperature  profile  of  the  atmosphere,  including  the  possible  existence  of  an 
inversion  layer.  Another  modification  takes  into  account  the  entrainment  of  dust 
particles  and  their  setting  at  their  size-dependent  terminal  velocity  relative  to  the 
gas  cloud. 

3B.  Wind-Controlled  Diffusion 

In  Figure  11  we  show  a  picture  taken  at  Ft.  Sill,  Ok.  during  the  Dust/Debris 
Tests  conducted  in  May  1978  by  Dugway  Proving  Ground  at  the  request  of  the  Project 
Manager  Smoke/Obscurant  Office,  Aberdeen  Proving  Ground,  MD.  Evident  in  the 
photograph  are  three  explosions  separated  in  time,  each  of  which  shows  a  strong  wind- 
shear  effect,  i,  e. ,  dust  near  the  ground  is  not  subjected  to  as  much  wind  advection  as 
dust  higher  in  altitude.  This  is  due  to  the  fact  that  the  wind  has  an  increasing  velocity 
with  altitude.  Also  evident  in  this  picture  at  the  leading  edge  of  the  dust  clouds  is  a 
puff  at  the  top.  This  may  be  associated  with  dust  that  was  entrained  in  the  initial 
buoyantly  rising  hot- gas  fireball. 

It  is  assumed  that  the  process  of  wind-controlled  diffusion  is  controlled  by  the 
differential  equation  of  diffusion,  the  action  of  turbulence  being  represented  on  a 
long-time  average  basis  by  the  gradient  transfer  relation.  The  diffusivities  in  the 
three  coordinate  directions  are  prescribed  consistent  with  the  stability  characteristics 
of  the  atmospheric  boundary  layer  (Ref.  12).  Aerosols  are  subjected  to  gravitational 
settling  (at  their  size-dependent  terminal  velocity)  in  addition  to  diffusion. 

The  initial  cloud  is  divided  into  horizontal  slices  of  circular  cross  section. 

Since  the  equations  are  linear  in  the  particle  concentration,  the  solutions  for  each 
slice  will  be  superimposed  to  produce  the  final  solution.  The  diffusion  of  the  gas  and 


Figure  11.  Example  of  wind-controlled  diffusion  of  dust  clouds  produced  by 
artillery  explosions  occurring  in  the  Dust/Debris  Test  conducted 
by  Dugway  Proving  Ground  at  Ft.  Sill,  Ok,  May  1978.  The  elon¬ 
gation  of  the  cloud  along  the  ground  is  due  to  the  effect  of  strong 
vertical  wind  shear. 

(Photograph  courtesy  of  Col.  Henry  R.  Shelton,  Mr.  Anthony 
Van  de  Wal,  and  Mr.  Eugene  Bowman,  PM  Smoke/Obscu rants 
Office,  Aberdeen  Proving  Ground,  Md.  ) 
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aerosols  contained  in  each  slice  is  then  described  by  solutions  of  the  diffusion  equation 
obtained  by  transform  methods.  Specifically,  we  consider  moments  of  the  particle 
concentration  function  C(x,y,  z,t),  viz.  , 


oo 

■  Jf 


C  dxdy 


where  x  is  the  coordinate  in  the  windward  direction.  The  diffusion  equation  may  be 


written  as 


§7  +  v  *  r  =  0 


Here  /^denotes  the  particle  flux  vector  with  components 

r  =  uc  -  d  ~ 

x  x  ax 


r  =  -  d  ~ 

y  y  ay 


r  =  -  v  C  -  D  (3c) 

z  g  z  az  ' 

Here  u  is  the  wind  velocity;  D  ,  D,,  and  D  are  the  diffusivities  along  the  x,  y,  and  z 

x  y  z 

directions;  v^  is  the  particle  terminal  settling  velocity.  It  is  important  to  note  that 

the  quantities  u,  D  ,  D  ,  and  D  are  all  allowed  to  vary7  with  z. 
x  y  z 

Multiplying  Eq.  (2)  by  appropriate  powers  of  x  and  y  and  then  integrating  over 
all  x  and  y  yields  the  moment  equation 

(at  0Z  Dz  az  ~  vg  az)  ^"pq  “  Pu  ^p-l,q  +  p^"^^x^p-2,q 
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Solutions  for  the  first  several  moments,  respectively,  describe:  the  mass  distribu¬ 
tion  of  dust  in  the  cloud  as  a  function  of  height  Z  and  time  t  (p  =  q  =  0);  its  mean 
along-wind  (x)  position  (p  =  1,  q  =  0);  its  mean  cross- wind  (y)  position  (p  =  0,  q  =  1); 
its  along-wind  spread  (p  =  2,  q  =  0);  and  its  cross-wand  spread  (p  =  0,  q  -  2). 

With  variables  made  dimensionless  with  respect  to  a  length  h  equal  to  the 
height  of  release,  a  velocity  jpz(h)/hj  ,  and  a  time  £h2/Dz(h)j  ;  then,  the  equation 
for  the  zeroth  moment  Tqq,  subject  to  a  diffusivity  =  Dz(h)Zn,  initial  conditions 

t  =  0  .  T00=  6  (Z-l)  (5- a i 

h 

and  boundary  conditions 

z  =  0  ,  ZnOT00/9Z)  =  0  (5-b) 

Z  ~  ~  ’  T00  “*  0  ’  (5-C) 


can  be  solved  in  closed  form  if  v  =0  to  obtain 
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where  ^n«2)/(2-n)  ^enotes  the  modified  Bessel  function  of  the  first  kind  of  order 
[(n-l)/(2-n)] . 

If  v  4  0,  a  closed  form  expression  for  the  zeroth  moment  is  obtained  by 

B 

heuristic  matching  of  two  asymptotic  solutions  obtained  by  the  following  process.  The 
equation 
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where  a  =  v  h/Dv(h),  is  recast  in  terms  of  the  independent  variables, 
g  z. 

£  =  Z  +  cvt  -  1  ,  T  =  t 


to  obtain 
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In  turn,  Eq.  (8)  is  approximated  by 
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in  £  »  err  ,  and  by 
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in  £  »  cut  .  The  solution  to  Eq.  (9-a),  subject  to  the  initial  and  boundary  conditions 
(Eqs.  (5-a,  b,  c) ),  is  readily  transposed  from  Eq.  (6),  viz., 
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The  solution  to  Eq.  (9-b),  subject  to  the  initial  condition  Eq.  (5-a),  the  requirement 
of  boundedness  as  |  —  ±  m  ,  and  the  restriction  (1  -  ar)  >  0,  takes  the  classical 
Laplace  form 
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where  ¥ = 


Z  >  0  and  r  <  a 


ar)n  dj.  A  heuristic  composite  solution  to  Eq. 
may  then  be  written  as 


(8),  valid  in 
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T00({,r,a)  =  f<5,ar)T^(«,T)  +  [l-  f(« ,  ar)]^  («,  r)  .  (11  > 


Here,  f(| ,  a t)  denotes  a  suitably  selected  weight  function  (see  below)  satisfying  the 
requirements  f  —  0  in  £  »  ar,  f  — *  1  in  £  «  ar. 

A  judicious  choice  of  the  function  f(| ,  ar)  is  obtained  upon  a  study  of  the  exact 
solution  of  Eq.  (7)  for  the  special  case  n  =  1,  i.  e. ,  the  case  of  neutral  stratification, 
where  the  diffosivity  varies  linearly  with  the  height  Z.  In  that  case,  the  change  of 
dependent  variable 

T00(Z,t)-=  z“T00(Z,t) 

leads  to  the  equation 


with  solution,  subject  to  the  initial  and  boundary  conditions  (Eqs.  (5-a,  b,c)), 

T00(Z,t)  -  Z°T00<Z,t)  [-^i]  I  [^]  .  (13! 

h 

Whereas  the  solutions  given  by  Eqs.  (10-a,b)  are  also  known,  the  weight  function 
f(| ,  ar)  is  readily  determined  upon  comparison  of  Eqs.  (11)  and  (13).  If  this  result, 
strictly  obtained  for  n  =  1,  is  adopted  throughout  the  practical  range  of  values 
0  <  n  ^  4/3  associated  with  stable  to  unstable  stratifications,  the  closed  form  solution 
Eq.  (11)  for  the  zeroth  moment  becomes  completely  defined.  The  availability  of  this 
solution  substantively  simplifies  the  determination  of  the  higher  moments  along  the 
lines  presented  below. 

Returning  to  Eq.  (4),  we  note  that  the  same  differential  operator  appears  in  the 

equations  for  all  the  moments.  However,  while  the  equation  for  Tqq  is  homogeneous, 

those  for  the  higher  moments  include  known  forcing  functions  R  .  Whereas  the 

pq 

initial  conditions  are:  t  =  0,  Tqq  «  6(Z  -  1),  and  T  =0;  the  solution  Eq.  (11)  for 
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Tqq  provides  the  Green's  function  required  to  construct  solutions  for  the  higher 
moments  by  superposition.  Specifically,  if  we  let 


G(4>£  ,  T  ,  T  ,  Of )  5 


then,  the  solution  for  the  general  moment  becomes 

or  'f 

Tp> q({ ,  r)  ~  I  d^'  f  Rp  q(|f.  T' )  Gg,  4',  T,  T  ,  O')  dT*  ,  (15) 

or  1  ^0 

which  can  be  readily  evaluated  by  numerical  integration, 

A  matrix  of  solutions  given  by  Eq.  (15)  covering  the  ranges  of  the  parameters 
n  and  a  associated  with  Pasquill  stability  categories  A  through  F  and  dust  particle  sizes 
D  ^  100  have  been  evaluated.  The  results  in  tabular  form  are  included  in  the 
DIRTRAN-I  code  together  with  interpolation  routines  which  permit  determination  of  the 
solutions  for  intermediate  values  of  the  parameters.  By  this  process,  the  calculation 
of  the  dust  concentration  within  a  given  dust  size  bin  at  a  prescribed  space /location 
and  time  achieves  the  same  simplicity  as  the  evaluation  of  a  Gaussian, 

A  typical  sequence  of  particle  concent  ration- length  product  (relative  line-of- sight 
concentration  integrals,  i.  e. ,  Cf  values)  contours  as  computed  by  the  DIRTRAN  model 
is  shown  in  Figure  12  for  a  sp’  re  of  initial  radius  5m  (corresponding  roughly  to  a 
155  mm  HE  artillery  shell  expl  *  >ion)  centered  at  x  -  0  and  z  =  5m.  The  times  plotted 
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are  5,  10,  30,  and  50  sec  after  detonation  in  Figure  12  for  purely  wind- controlled 
case  and  in  Figure  13  for  the  initially  buoyancy  controlled  case.  In  both  cases,  a 
wind  speed  of  3  m/sec  (measure  at  height  of  10m  above  the  ground)  was  used. 

Figure  12  is  the  result  of  calculations  for  the  case  of  small  particles,  i.  e. ,  no 
gravitational  settling.  Figure  13  shows  the  results  for  particles  of  nominal  diameter 
50  /j m  and  a  nominal  gravitational  settling  velocity  of  18  cm /sec.  Figure  14  shows 
the  superposition  of  an  initially  buoyancy- dominated  case  (which  is  subsequently  wind 
controlled)  with  the  immediately  wind  dominated  case,  both  far  the  situation  of  small 
particles.  Noteworthy,  in  Figure  14  is  the  development  of  a  stretch ed-out  cloud 
hugging  the  ground,  quite  similar  to  that  manifested  in  the  data  of  Figure  11. 

In  conclusion,  we  submit  that  the  DIRTRAN  model  reproduces  the  salient  features 
of  dust  clouds  in  a  battlefield  scenario.  Detailed  validation  with  data  from  the  DIRT  I 
and  Graf  II  test  series  is  being  actively  pursued  at  the  present  time. 
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Wind _  Wind  Speed:  3m/sec  (measured  at  10m  height) 

Direction  ""  Wind  Power  Law  Exponential:  1/7 

Initial  Sphere  Centered  at  x^O,  5m  height) 


t*30  sec 


t  *  50  SEC 


Distance,  meters 


Figure  12.  Wind  controlled  aerosol  dispersion.  Time  sequence  of  contours 

of  line-of- sight  integrals  of  dust  concentration  per  unit  area  (i.  e. , 
Cf  values)  computed  by  the  DIRTRAN  model  for  the  case  of  no 
initial  buoyancy.  Note  the  effect  of  wind  shear,  particles  are 
assumed  small  (i.  e. ,  negligible  gravitational  settling).  The  three 
concentric  contours  represent  respectively,  0.1,  0.05,  and  0.025 
fractions  of  the  initial  average  line-of-sight  integral  of  particle 
concentration. 


DIRTRAN  CALCULATION 


Wind _ Jfiwo  Speed:  3m/sec  (measured  at  10m  height) 

Direction  *Wind  Power  Law  Exponential:  1/7 

Initial  Sphere  Centered  at  x=0,  5m  height) 


T*10  SEC 


t*50  sec 


T*5  sec 


Distance,  meters 

Figure  13.  Buoyancy  diffusion  of  aerosols  with  wind  and  gravitational 
settling  (50  jiin  diameter  particles).  Note  that  for  t  ^  30 
sec,  the  dust  cloud  has  sufficiently  diffused  that  the  0. 1 
fraction  contour  no  longer  exists. 


DIRTRAN  CALCULATION 


Figure  14.  Superposition  of  wind  and  buoyancy  controlled  diffusion  of 
small  particles.  Note  the  strong  wind  shear  effect  with 
altitude,  an  important  phenomenon  for  near-ground  line- 
of-sight  obscuration  effects. 
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4.  IR  CLUTTER 


In  addition  to  total  obscuration  of  a  target  of  interest  by  battle  field- induced  con¬ 
taminants,  there  is  the  serious  problem  of  battlefield  clutter  produced  by  these  par¬ 
ticulates  as  well  as  associated  turbulence.  In  this  situation,  there  may  be  so  much 
spatial  and  temporal  variation  in  a  detector  field  of  view  due  to  Cj?  (aerosol 
concentration- length  product)  or  thermal  changes  that  a  target  may  not  be  readily 
discernible.  This  can  happen  even  though  the  target  may  not  be  completely  obscured. 
This  clutter  can  result  in  visibility  holes  at  infrared  wavelengths  and  scintillation  at 
mm  wavelengths.  This  was  especially  evident  in  the  Graf  II  tests  with  tank  night 
sights  operating  in  the  8-12  /^m  region.  At  mm  wavelengths,  results  from  the 
DIRT  I  tests '  showed  that  scintillation  could  be  a  factor  of  three  greater  than  the 
average  attenuation.  These  marginal  seeing  conditions  pose  serious  threats  to 
electro-optical  systems,  and  an  understanding  of  the  reasons  for  and  the  probability 
of  their  occurrence  is  essential  if  one  is  to  ascertain  system  effectiveness.  We  sug¬ 
gest  here  that  the  existence  of  battlefield  clutter  may  be  attributed  to  the  interplay  of 
barrage  intensity  (be  it  high  explosive  or  smoke  rounds)  and  large  atmospheric  eddies 
of  characteristic  period  and  size. 

Inspection  of  available  transmissometer  data  forcibly  exhibits  the  existence  of 
distinct  transmission  windows  (visibility  holes)  in  battlefield-induced  obscurants  under 
marginal  seeing  conditions.  The  visibility  holes  under  these  barrage  conditions  can 
be  explained  by  the  presence  of  large  atmospheric  turbulent  eddies.  In  the  temporal 
domain,  these  transmission  events  have  typical  individual  durations  on  the  order  of 
ten  seconds  and  repetition  rates  on  the  order  of  one  minute.  In  the  spatial  domain, 
the  characteristic  wavelength  is  of  the  order  of  150m  between  events,  and  the  scale 
of  each  event  is  on  the  order  of  30m.  The  data  in  Figure  15  show  the  visibility  holes 
observed  during  the  Ft.  Knox  TAPATS  test  conducted  by  NV&EOL  in  July  1978.  The 
optical  transmission  level  between  these  holes  is  relatively  monotone  in  time;  the  ac 
component  of  the  transmission,  however,  is  dominated  by  large  atmospheric  eddies 
which  are  of  well  defined  scale  (size)  and  frequency  and  not  by  a  spectrum  of  random 
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Transmission  (l) 


VISIBILITY  HOLES-  INJiVSlQLJAFAIl  BARREL  --DAYTIME,- 26  iOLPOZB 

(Ft.  Knox,  Ky.)  Wavelength  Band:  8  -  12  n* 

Artillery:  155  mm  HE 
Barrage:  0,7  Rounds/Km/sec 
Impact  Area:  100  m  x  300  m 
Atmosphere:  Unstable 

WlNDSPEED’  4.5  KNOTS 


Figure  15.  Effect  of  large-eddy  turbulence  on  visibility  holes,  periodicity 
of  42  sec  and  duration  of  10  sec  per  event  is  evident.  This  can¬ 
not  be  explained  by  the  barrage  itself  which  was  an  average  cf 
one  round  every  five  seconds  in  rapid  random  fire  from  an  artil¬ 
lery  battery  of  six  guns.  (The  probability  is  negligible  that  the 
random  rounds  would  land  in  the  same  place  with  respect  to  the 
line  of  sight  at  precisely  the  same  interval  and  duration  at  three 
different  times  during  the  barrage.  In  addition,  there  is  another 
hole  (an  obscuration  hole)  occurring  after  the  barrage  and  which 
is  of  the  same  periodicity  and  duration  as  the  other  three  holes. ) 
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(turbulent)  scales  and  frequencies.  A  deterministic  model  descriptive  of  the  nature, 
scale,  and  frequency  of  these  dominant  eddies  is  required  to  predict  the  occurrence 
of  the  windows  and,  therefore,  the  nature  of  clutter  to  be  encountered  in  a  battlefield 
under  the  conditions  of  marginal  viewing. 

The  influence  of  large  eddies  on  specific  measurements  of  obscuration  is  also 
clearly  manifested  in  other  currently  available  data.  The  data  in  Figures  16  to  19 
are  from  the  DIRT  I  tests.  The  pictures  on  the  fire  smoke  obscurant  clouds  show 
markedly  different  cloud  patterns  with  different  lengths  and  frequencies  as  mentioned 
above.  A  comparison  of  the  pictures  also  suggests  that  validation  of  any  obscuration 
model  will  be  prevented  by  the  limited  number  of  data  available  and  the  consequent 
inability  to  determine  the  Mattendant  statistics",  we  wish  to  point  out  that  the  narrow 
bandwidth  of  the  large-eddy  effects  (i.  e. ,  the  presence  of  a  well  defined  characteristic 
frequency  for  a  given  meteorological  condition)  clearly  shows  two  options  for  model 
validation:  1)  the  repetition  of  the  same  measurement  for  a  statistically  meaningful 
number  of  events  (>  10  to  100)  and  subsequent  ensemble  averaging  of  the  results  at 
fixed  meteorological  condition  and  obscurant  source;  or  2)  the  analytical  prediction 
of  the  variance  in  the  observations  produced  by  the  large  turbulent  eddies. 

We  conclude  that  an  understanding  of  the  interplay  between  barrage  intensity  and 
large  atmospheric  eddies  of  characteristic  period  and  size  will  provide  insights  to 
the  problem  of  battlefield  clutter,  scintillation,  and  the  existence  of  visibility  holes. 

At  present,  we  are  directing  efforts  toward  resolving  this  problem. 


Figure  16.  Example  of  visibility  holes  produced  by  the  interplay  of  large-edd 
turbulence  effects  and  fire  smoke  generation  on  obscuration. 
(Tests  performed  during  the  DIRT  I  series.  ) 


irt 


Firt  :  lufsi  J  oil;  rubber;  pn 
ST  A I  I  'A :  F-  r  -  *>  TIX1 


>** 


*sbb 

■&x£. 

lilt 

v 

<&- 1 

Figure  17. 


Two  pictures  of  buoyantly  rising  smoke  separated  in  time  by 
approximately  40  sec.  In  the  lower  event,  the  large- eddy 
turbulence  is  sufficient  to  completely  offset  (nullify)  the  rise 
of  the  smoke  due  to  buoyancy  from  the  hot  gases  produced  in 
the  fire.  In  this  situation,  a  visibility  hole  is  created  at  an 
altitude  of  about  70m.  (Data  from  DIRT  1. 1 
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Figure  18  Example  of  biftircation  of  DIRT  I  fire  smoke  clouds  due  to  large 
atmospheric  eddies  separated  in  time  by  approximately  20  sec. 
A  visibility  hole  between  the  two  plumes  at  an  altitude  of  40m  is 
evident.  Obscuration  effects  at  ground  level  are  also  evident. 
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Figure  19.  View  of  approximately  the  same  MwishboneM  fire  smoke  event 
(temporal  separation  of  less  than  2.  5  sec)  from  two  different 
aspects  (angular  separation  of  64.  2°).  In  the  upper  picture 
the  large-eddy  turbulence  has  more  than  offset  the  buoyant  rise 
of  the  cloud,  i.  e. ,  the  smoke  has  been  forced  back  to  the  ground 
creating  visibility  blockage  along  the  ground  (zero  altitude)  at  a 
distance  of  40m  from  the  fire  smoke  source.  A  visibility  window 
is  evident  20m  laterally  from  the  source  near  the  center  of  the 
"wishbone"  smoke  distribution.  In  the  lower  picture,  no  visi¬ 
bility  windows  are  evident  whatsoever  due  to  the  different 
aspect  angle.  (Data  from  DIRT  I. ) 
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5.  SUMMARY  AND  CONCLUSIONS 


Fbr  the  dust  aerosol  source  term,  we  conclude  that  dust  particles  inside  and 
outside  of  the  explosion- produced  fireball  have  different  diffusion  mechanisms  and 
obscure  in  different  regions  of  the  atmosphere  as  evidenced  by  DIRT  I  data;  they, 
thus,  need  to  be  accounted  for  separately  and  then  superimposed.  We  also  conclude 
that  the  large  differences  in  obscuration  between  the  Ft.  Knox  TAPATS  (Threat 
Artillery  Preparation  Against  Tank  Sights)  test  and  the  Grafenwohr  Graf  II  test  may  be 
attributed  to  the  existence  of  thick  (^25  cm)  vegetative  cover  (sod)  at  Grafenwohr, 
resulting  in  a  smaller  amount  of  entrained  soil  particles.  For  the  atmospheric  dif¬ 
fusion,  we  discuss  provisions  in  the  DIRTRAN  model  for  analytic  solutions  (which  run 
fast  on  the  computer)  to  the  atmospheric  diffusion  equations  which  account  for: 
vertical  wind  shear  (important  for  near-ground  linc-of-sight  obscuration  effects); 
variable  vertical  diffusivities  with  altitude  (important  for  different  atmospheric  stabil¬ 
ity  categories);  and  gravitational  settling  (important  for  large  dust  particles).  It  is 
observed  that  these  solutions  realistically  model  the  time-dependent  cloud  shape  of 
battlefield  aerosols.  Finally,  we  show  that  the  existence  of  battlefield  clutter  which 
occurs  in  a  barrage  is  the  result  of  the  interplay  of  barrage  intensity  and  large  atmos¬ 
pheric  eddies  of  characteristic  period  and  size.  This  has  a  strong  dependence  on 
stability  category,  and  produces  visibility  holes  at  IR  wavelengths  and  scintillation  at 
mm  wavelengths. 
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